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We propose a model for fiigfi-Tc superconductors, valid for < 5 < Ssc, tfiat includes both the spin 
fluctuations of the Cu^"'' magnetic ions and of the 0~~ doped holes. Spin-charge separation is taken 
into account with the charge of the doped holes being associated to quantum skyrmion excitations 
(holons) of the Cu++ spin background. The holon effective interaction potential is evaluated as a 
function of doping, indicating that Cooper pair formation is determined by the competition between 
the spin fluctuations of the Cu^''' background and of spins of the O doped holes (spinous). The 
superconducting transition occurs when the spinon fluctuations dominate, thereby reversing the sign 
of the interaction. At this point {5 — 5sc), the theory is supersymmetric at short distances and, 
as a consequence, the leading order results are not modified by radiative corrections. The critical 
doping parameter for the onset of superconductivity at T = is obtained and found to be a universal 
constant determined by the shape of the Fermi surface. Our theoretical values for Ssc are in good 
agreement with the experiment for both LSCO and YBCO. 



PACS number(s): 74.72.Bk, 74.25.Ha 

Introduction. High-temperature superconducting 
cuprates have a very rich and complex phase diagram 
whose understanding is an important issue. In the un- 
derdoped region, for instance, a wide variety of physi- 
cal phenomena like Neel and metal- insulator transitions, 
transport (non-Fermi liquid) anomalies, the occurrence 
of a spin-pseudogap, absence of a sharp quasiparticle 
peak (spin-charge separation), etc., have inspired a large 
amount of theoretical and experimental work for about 
fifteen years In spite of that, even the nature of the 
ground state and of its elementary excitations has not 
yet been fully determined and many different pictures 
are availabe, ranging from a dimerized ground state with 
spin-Peierls or valence-bond order |^ until the so called 
staggered- flux (d-wave) phase 

Another fundamental point yet to be understood is 
the mechanism of charge pairing leading to the super- 
conducting transition. In connection to this, we stress 
that it is by now well established that antiferromagnetic 
spin correlations play an important role in the dynamics 
of the system, even after the destruction of the Neel state. 
Indeed, different spin-fluctuation models have been suc- 
cessfully used to explain the observed spectral weight in 
ARPES data in the quantum disordered phase of high-Tc 
materials [Q, as well as other anomalies [||. Moreover, 
the idea of spin-fluctuation induced charge pairing and 
superconductivity has been used recurrently 

In this work we propose a model for high-Tc cuprates 
valid for dopant concentrations ranging from zero up to 
the superconducting transition, < 5 < 5sc^ that takes 
into account the spin fluctuations of the Cu"'""'" magnetic 
ions and of the 0~~ doped holes on different footing, as 
suggested by the different temperature dependences of 



the NMR relaxation rates for the Oxygen and Copper 
spins 0. Our model also incorporates spin-charge sep- 
aration as follows. The charge of the dopants intro- 
duced in the Cu02 planes is associated to skyrmion quan- 
tum spin excitations of the Cu"'""'" background (holons) 
which, in the Neel phase appear as finite energy defects 
closely related to their classic counterparts, whereas in 
the quantum disordered phase are nontrivial zero energy 
purely quantum mechanical excitations. The spin of the 
doped holes, on the other hand, is represented by charge- 
less, massless Dirac fermion fields (spinons) Q. We then 
calculate the effective interaction potential between the 
quantum skyrmion topological excitations, as a function 
of doping, in order to study charge pairing. It becomes 
clear that Cooper pair formation at zero temperature is 
controlled by the competition between two different con- 
tributions to the quantum skyrmion effective interaction 
energy; one coming from the spin fluctuations of Cu"'""'" 
magnetic ions and the other from the corresponding fluc- 
tuations of the spins of the doped holes (spinons). The 
superconducting transition occurs when the latter domi- 
nates, thereby reversing the sign of the effective quantum 
skyrmion (dopant charge) interaction potential, at short 
distances, from repulsive to attractive. Interestingly, the 
model becomes supersymmetric at short distances, pre- 
cisely at the onset of superconductivity, thereby making 
our one-loop results robust against radiative corrections, 
in accordance to Witten's theorem The critical 

value of the doping parameter for the superconducting 
transition at zero temperature, Ssci is universally deter- 
mined by the shape of the Fermi surface and is in good 
agreement with experiment for both La2-a:Sra:Cu04 and 
YBa2Cu306+a; compounds. The model also correctly 
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predicts the zero temperature magnetization curves as a 
function of the critical doping in the AF ordered phase. 

The model. In previous works | [Tl| ], we have proposed 
a model for doping quantum Heisenberg antiferromag- 
nets that successfully described the magnetization curves 
and the AF part of the phase diagrams of the two best 
studied high-Tc compounds, namely LSCO and YBCO. 
One of the important consequences of that model is the 
observation that each hole added to the Cu02 planes 
creates a skyrmion topological defect, as has been pro- 
posed earlier The dopant charge, in particular, is 
attached to the skyrmion charge and its dynamics be- 
comes totally determined by the quantum skyrmion cor- 
relation functions. In the Neel ordered Mott insulating 
phase, the skyrmions have a finite excitation energy and 
this reflects the existence of a gap for charge that can 
be associated to the breakdown of translational invari- 
ance of the lattice. The model proposed in however, 
is restricted to the antiferromagnetic part of the phase 
diagram, where S < Saf- Nevertheless, we shall pur- 
sue the picture in which skyrmions are in general the 
charge carriers of the doped holes. In particular, we 
shall exploit this idea in the quantum disordered phase, 
^AF < (5 < 5scy where the skyrmions are purely quan- 
tum mechanical and have zero energy. This is again con- 
sistent with charge response experiments, like in-plane 
optical conductivity which indicates the absence of 
a gap for charge excitations in the quantum disordered 
metallic phase. 

Let us consider the zero temperature Euclidean parti- 
tion function 
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and we use units in which Ti = cq = 1. In the above 
expression zj,Zi, i — 1,2, are Schwinger boson fields re- 
lated to the local spin density of Cu''"*" ions through S — 
z\(TijZj, and ipl^ipa, a =t, i, are chargeless 2-component 
Dirac spinor fields (spinous) describing the local spin den- 
sity of the doped O holes through Sh = iplcTabipb- As 
usual, Ai_i is the Hubbard-Stratonovitch field and go the 
bare coupling constant of the CP^ model. The constant 
q measures the strength of the coupling of spinous. 

Spin-charge separation is manifested in our model 
through the fact that the massless Dirac fermions carry 
the spin of the doped holes, whereas all information 
about their charge is carried by the quantum skyrmion 
excitations (holons) created out of the Schwinger boson 
background [Q. The description of spinous as massless 
Dirac fields arises naturally in the continuous limit of mi- 
croscopic models B . The full treatment of the quantum 



skyrmions of the theory described by (^, on the other 
hand, has been carried out in [l^ . 

Before we proceed, it is important to determine how 
the doping dependence will be introduced in our theory. 
As we explained above, we may identify in principle, at 
least two phases in the model given by (^), at T = 0. An 
ordered Neel phase, go < gc, for which there is a nonzero 
spin stiffness ps = 1/ go — 1/ gc > {gc is the quan- 
tum critical coupling) and a quantum disordered phase, 
go > gc, in which the Schwinger bosons Zi acquire a mass 

oc [1/gc — l/ffo]i and ps = 0. As we explained in 
the whole dynamics of the in-plane dopant charge is iden- 
tical the quantum dynamics of skyrmions. The quantum 
skyrmion correlation function corresponding to (|^) has 
been evaluated in (l4| , in the ordered phase, giving 



{p{x)p\y)) = 
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where /i^ is the quantum skyrmion creation operator. 
Conversely, for the theory studied in ||T^ the correspond- 
ing correlator was found to be 
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where the expressions for Ps{5) and ol{5) have been de- 
termined in |11|. In particular. 
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with n = 1 for YBCO and n = 4 for LSCO, the factor of 
four being a consequence of the existence of four branches 
in the Fermi surface for this compound, see discussion in 
The ps{5) fimction is given by p{5) = p(0)[l - A5\ 
for YBCO and p{5) = p(0)[l - B5 - CS^Y'"^, for LSCO, 
and again the different behavior being ascribed to the 
form of the Fermi surface in each case [0. The con- 
stants A,B and C have been evaluated from first prin- 
ciples in jllj. In (^), uem is the electromagnetic fine 
structure constant and accounts for the contribution of 
the electromagetic interaction of the doped holes to the 
skyrmion correlation function. Examining we see that 
actually this term can be neglected when compared to 
the first one. In order to obtain the (5-dependence of 
the spin stiffness ps and of the spinon coupling q in our 
model (^), we now match the two correlation functions 
in @ and (0) (ordered phase), obtaining ps = p{5) and 
q = [ ^l^^^g ]~^(n^), where we have already neglected the 
electromagnetic part. We immediately conclude that the 
sublattice magnetization in the ordered phase is given by 
M(5) = ^/p(J). From this we can readily obtain 5af 
from p(5af) — 0, see also [0. For 5 > Saf, on the other 
hand, where = 0, we assume that the expression for 
q{S) still holds. 

Cooper pair formation. Let us now investigate the con- 
ditions for skyrmion pairing and consequent formation of 



2 



Cooper pairs, by analyzing the effective interaction po- 
tential between quantum skyrmions in the quantum dis- 
ordered underdoped phase. For this purpose, we intro- 
duce the skyrmion current = ^^^'^"^Sc^^ through 
the identity 
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Integrating over z|, Zi and ^p^jij^a, we obtain, at one- loop 
level, the effective action 

(7) 

where the kernel S(a; — y) has Fourier transform given by 
Y.{p) = nB(p) + n^b) where 



and 



1 

2^ 



m 



2p 



-arctan 



— ^arctan - — 

V TO / p'^ \ 2ni / 



Uf{p) = 



27r 



(8) 



(9) 



These two terms are, respectively, the contributions to 
the vacuum polarization coming from the complex scalar 
fields Zi (Schwinger bosons) and fermions ip (spinous). 
In (H), TO is the mass of the z^-fields (spin-gap) in the 
quantum disordered phase, where 6 > Saf (ffo > 5c) 

In order to obtain the effective skyrmion action, we 
use an exponential representation for the (5-function in (|^) 
and integrate over the corresponding Lagrange multiplier 
field and Afj,. The result is 
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where the Fourier transform of the kernel is given by 
S(P) 
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From expression (|T^ we can readily obtain the effec- 
tive interaction energy between static skyrmions. This is 
given, in real time, by 

Hi = -^i^TTf I d^x I d2yp(x)/^(x - y)p(y), (12) 

where /o(x) = j7o(x, 0) is the dopant charge density and 



A'(x-y) = J d2p/(27r)2 e'P'^'^-y) S(p,0). 

For two charges at positions xi and X2, we have 
p(x) = (5(^) (x — Xi) -t- S^"^^ (x — X2). After discarding self- 
interactions, we obtain the effective interaction potential 
for static charges, namely 



V^(xi-X2)= yd2pl](p,0) e'P (^^-^^) 
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It is well known that in high-Tc superconducting 
cuprates. Cooper pairs form in such a way that the two 
charges are localized in space. Only the short distance 
behavior of the interaction potential, therefore, is rele- 
vant for Cooper pair formation. In this limit (large |p|) 
we have 



F(xi -X2) ^ Jd^p 
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The first contribution inside the square brackets in the 
above expression corresponds to the Cu"*""*" spin fluctu- 
ations (Schwinger bosons) while the second corresponds 
to fluctuations from O^^ spins (spinous). We see that 
for small doping, < 1, the potential is allways repul- 
sive and there is no charge pairing. For > 1, on the 
other hand, the interaction potential becomes attractive 
and charge (skyrmion) pairing occurs. Consequently, we 
conclude that the critical doping for the onset of super- 
conductivity is determined by the condition 
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Let us remark that if we had considered a system of 
spinous solely, without including the Cu++ background, 
we would have obtained that the interaction potential 
(|l5| ) would always be attractive for any g 7^ 0, at zero 
temperature, and 6sc = 0- This is what happens in the 
mean field phase diagram of Kotliar and Liu for the t — J 
model [[|. We see that the primer effect of considering 
the Cu"*""*" spin background is to shift the value of Ssc to 
the right of the phase diagram, which is actually what is 
observed experimentally. 

Comparison with experiment. From the expression of 
q in terms of (5 (g = 2^/a{S), see (|^)), we may infer 
that 6sc is an universal constant, only determined by 
the shape of the Fermi surface. We see, in particular, 
that 5%c^^ = '^^scF^i a relation that is verified by ex- 
periments, if we take in account the relation between i5 
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and the stoichiometric doping parameter x, namely 5 = x 
for LSCO and 5 = a;-0.20 for YBCO. Another prediction 
of our model is that compounds with similar Fermi sur- 
faces should have the same superconducting critical dop- 
ing 5sc- From (|l|) and (|), we calculate S^S^° = 0.318 
and Sgff"^ ~ 0.079, which have a fairly good agreement 
with experiment. We show below that taking in account 
the presence of disorder, we can obtain better values for 
these critical doping parameters. 

Disorder. Disorder may be modelled in the ordered 
Neel phase of a doped antiferromagnet by considering a 
continuous random distribution of spin stiffnesses . If 
we introduce a Gaussianxjp — Ps\'^~^ distribution, with 
exponentially suppressed magnetic dilution, in the orig- 
inal model |]ll| used to describe the antiferromagnetic 
phase, we obtain a correction for (^, namely a{5) 
a'{5) = a{6) -I- v ||l^. Choosing u — ^ ioi both com- 
pounds, we get 



LSCO and YBCO compounds. We stress that the pairing 
must be between skyrmions and not between skyrmion 
and anti-skyrmion, so that the total electric charge of the 
pair is nonzero. 

At finite temperatures, the pairing shall no longer oc- 
cur for q^{5) ~ 1. There will be finite temperature cor- 
rections for both Hb and Tip and, since supersymmetry 
is broken at any finite temperature, the fermionic and 
bosonic contributions for the interaction potential should 
cancel a.i q ^ 1 {5 ^ 5sc)^ in agreement with experimen- 
tal results. We are presently investigating this point. 

This work has been supported in part by FAPERJ and 
PRONEX - 66.2002/1998-9. E. C. M. was partially sup- 
ported by CNPq and M. B. S. N. by FAPERJ. 
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Now, the critical doping parameters at T = become 
(Jl'sco ^ 0.225 and 5^$^° = 0.056, corresponding to 
x^^C'O ^ 0.425 and x^^^° = 0.056, which are in good 
agreement with experiment. Notice that with a single 
choice for the disorder distribution we correctly obtained 
the critical dopings for both YBCO and LSCO. 

Supersymmetry. Let us observe now a remarkable 
fact. At short distances, when the mass of the Schwinger 
bosons, m, may be neglected, our model becomes su- 
persymmetric precisely at the point q{5sc) = 1, where 
the superconducting transition occurs. Supersymmetry 
relates Schwinger bosons Zi and spinous ij^a and that is 
why the contributions of both to the holon (skyrmion) in- 
teraction potential are identical but with opposite signs 
at this point. An important consequence is that our one- 
loop derivation of the holon effective interaction potential 
and critical dopings are unchanged by radiative correc- 
tions. Indeed, we have actually checked that the con- 
tribution of these corrections to the short distance be- 
havior of the effective holon interaction potential ( p^ ) is 
subdominant and can be neglected. This can be under- 
stood on general grounds, as a result of Witten's theorem 
, which states that supersymmetry cannot be broken 
perturbatively. 

Conclusions. We have calculated the effective interac- 
tion potential between holons in a spin-charge separated, 
spin fluctuation model for hith-Tc cuprates. We have 
shown that Cooper pair formation and superconductiv- 
ity is determined by the competition between the spin 
fluctuations of the CU++ antiferromagnetic background 
and the spin fluctuations of the doped O holes, in the 
undcrdoped regime. Our prediction of the critical doping 
for the onset of superconductivity at zero temperature, 
6sc, is in good agreement with experiment for either 



[2] 
[3] 

[4] 
[5] 
[6] 



[7] 
[8] 

[9] 



[10] 
[11] 



[12] 



[13] 
[14] 
[15] 
[16] 



P. W. Anderson, The Theory of Superconductivity in the 
High-Tc Cuprate Superconductors, Princeton University 
Press, (1997); W. Brenig, Phys. Rep. 251, 153 (1995). 
P. W. Anderson, Science 235, 1196 (1987); N. Read and 
S. Sachdev, Phys. Rev. B 42, 4568 (1990). 

D. Poilblanc, E. Dagotto and J. Riera, Phys. Rev. B 43, 
7899 (1991); D. A. Ivanov, P. A. Lee and X-G. Wen, 
Phys. Rev. Lett. 84, 3958 (2000); P. W. Leung, Phys. 
Rev. B 62, R6112 (2000). 

W. Ranter and X-G. Wen, Phys. Rev. Lett. 83, 3871 
(2001); M. Franz and Z. Tesanovic, Phys. Rev. B 63, 
064516 (2001). 

A. V. Chubukov and J. Schmalian, Phys. Rev. B 57, 
R11085 (1998); J. Schmalian, D. Pines and B. P. Sto- 
jkovic, Phys. Rev. Lett. 80, 3839 (1998). 
T. Moriya, K. Ueda and Y. Takahashi, J. Proc. Soc. 
Jpn. 49, 2905 (1990); P. Monthoux, A. Balatsky and D. 
Pines, Phys. Rev. Lett. 67, 3448 (1991); T. S. Nunner, 
J. Schmalian and K. H. Bennemann, Phys. Rev. B 59, 
8859 (1999). 

M. Takigawa et al, Phys. Rev. B 43, 247 (1991). 

G. Baskaran and P. W. Anderson, Phys. Rev. B 37, 580 

(1988) . 

J. B. Marston and I. Affleck, Phys. Rev. B 39, 11538 

(1989) ; P. A. Lee, N. Nagaosa, T-K. Ng and X-G. Wen, 
Phys. Rev. B 57, 6003 (1998); D. H. Kim and P. A. Lee, 
Ann. Phys. 272, 130 (1999). 

E. Witten, Nucl. Phys. B202, 253 (1982). 

E. C. Marino, Phys. Lett. A 263, 446 (1999); E. C. 
Marino and M. B. Silva Neto, Phys. Rev. B (2001), in 
production (3ond-mat/0008186). 



P. B. Wiegmann, Phys. Rev. Lett. 60, 821 (1988); B. 
I. Schraiman and E. Siggia, Phys. Rev. Lett. 61, 467 
(1990); R. J. Gooding, Phys. Rev. Lett. 66, 2266 (1991); 
S. Haas, et al, Phys. Rev. Lett. 77, 3021 (1996). 
T. Timusk, et al, Physica C 162, 841 (1989). 
E. C. Marino, Phys. Rev. B 61, 1588 (2000). 
G. Kotliar and J. Liu, Phys. Rev. B 38, 5 142 (1988). 
E. C. Marino, preprint :ond-mat/0106514. 



4 



